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Abstract 

We investigate the sensitivity to new physics of the process e+e“ —>■ tt when the top polar¬ 
ization is analyzed using leptonic final states e+e“ —)■ 11 —>■ Rl~b b vivi. We first show that the 
kinematical reconstruction of the complete kinematics is experimentally tractable for this pro¬ 
cess. Then we apply the matrix element method to study the sensitivity to the Vti coupling {V 
being a vector gauge boson), at the tree level and in the narrow width approximation. Assuming 
the ILC baseline configuration, '/S = 500 GeV, and a luminosity of 500 fb“^, we conclude that 
this optimal analysis allows to determine simultaneously the ten form factors that parameterize 
the Vti coupling, below the percent level. We also discuss the effects of the next leading order 
(NLO) electroweak corrections using the GRACE program with polarized beams. It is found 
that the NLO corrections to different beam polarization lead to significantly different patterns 
of contributions. 


‘Proceeding for the series of the TYL/FJPPL workshops on Top Physics at the ILC held at KEK (November 
2013), LPNHE Paris (March 2014), and LAL Orsay (March 2015). 
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1 Introduction 


Top quark was discovered in 1995 as the 6 th quark. Its heavy mass, which is orders of magnitude 
higher than those of the other quarks, was predicted a decade before by indirectly measuring the 
radiative processes which receive a contribution from top quark propagating in the loop, which 
increases as the mass becomes higher. Today, a more precise measurement of the top quark mass is 
desired in order to discover a small deviation from the electroweak theory indirectly through loop 
corrections. Furthermore, top quark mass being close to the electroweak scale, one can speculate 
that the top quark may play a special role for the electroweak symmetry breaking, namely in the 
new physics models. 

On the other hand, the top quark has another peculiarity: its decay time is so short that it does not 
hadronize. This leads to a great advantage: the top quark is the only quark whose spin polarization 
can be studied. In that sense it is akin to the tau, whose production and polarization was studied 
at LEP and SLC, but with two remarkable qualitative differences due to its large mass: all four 
intermediate polarized states: ti,iji , , tjitji , tjiii are produced, and interfere, therefore 

creating a rich pattern bearing the imprint of the top mass; the VF’s are produced on-shell, and 
their decay can be fully reconstructed. The study of polarization allows us to probe the chirality of 
the interactions between the top quark and the gauge bosons 7 and Z or any new particles beyond 
the SM. Furthermore, the top polarization measurement also provides us an opportunity to study 
the CP violating interactions, in a particularly clean manner, unhampered by hadronic effects. The 
use of the top polarization can be done at e"*“e“ colliders, such as ILC [ 2 ], where massive numbers 
of top anti-top pairs can be produced. It is well known that lepton colliders, without any hadronic 
effect in the initial state, can provide a clean environment to study the top quark properties [2]. In 
particular, the option of the ILC proposal to use polarized beams is quite interesting as the initial 
state beam polarization and final state polarization exhibit strong correlation, which can be used 
as an additional information to study the top quark interactions. 

The top quark polarization pattern can be reached by measuring the angular distribution of the 
top quark decay products. For top and anti-top quarks decaying into the leptonic final states 
e'’“e“ —)• t t —>■ l~^l~b b viDi = e~^e~,e~^fj,~, fj,~) it is known that the full kinematics 

(hence including the top directions) can be obtained from the lepton and bottom momenta [5l[7]. 
In section we first demonstrate that, working at tree level, one can determine the so-called form- 
factor of the tt production with a precision down to the few per-mill level. In this analysis, we 
use the matrix element method which is the most efficient method when all the kinematics can be 
reconstructed. 

The obtained result is quite interesting. Thus, it appears feasible to investigate further beyond the 
tree level, namely including the 0(a) and the 0{as) corrections. It has been a long time since a 
large electroweak NLO correction for e'’'e“ —>• ti was recognized [3] and confirmed independently 
in [Ij. However, this effect has not yet been taken into account for the top polarization study. Since 
the foreseen experimental errors can be as small as per-mill level, the theoretical error coming from 
the higher order corrections has to be evaluated very carefully. 
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The remaining of the sections is organized as follows. In section we first revisit the feasibility 
of the kinematical reconstruction of the fully leptonic final states. In section we perform a 
sensitivity study in terms of the effective form factors which describe the ttZ/'j coupling. Section 
[^includes the discussion on the electroweak NLO corrections, and the conclusions and the future 
prospects are given in section 

In the following we assume a luminosity of 500 fb~^ collected at a fixed center of mass energy of 500 
GeV with an equal sharing in luminosity between the two beam polarizations of the ILC, which 
corresponds to about 10^ events with = —80% and Vf,+ = +30%, and 6 x 10^ events with 
Vf,- = +80% and 7^g+ = —30%. No acceptance cuts are applied; the reconstruction is assumed to 
be perfect; and detector and physics effects (ISR, hadronization, backgrounds, etc.) are ignored. 


2 Kinematical reconstruction of the fully leptonic final state: test 
with LO Monte Carlo events 




Figure 1: The plane of the top angles in the laboratory frame; </>( (horizontal axis), 9t (vertical axis). The 
two close curves represent the solutions of the kinematical equations for a particular event, randomly chosen. 
On the left panel, the two countours are obtained by imposing that the computed energy of the lepton in 
the rest frame of the W is equal to tovv/ 2. On this panel, the ellipsoidal contour on the left (resp. right) 
corresponds to I (resp. 1). These two countours cross in two points, which are the two solutions of the 
kinematical reconstruction. On the right panel, the two small additional ellipsoidal contours are drawn by 
using the +jet energy information: they correspond to the domains allowed at 99%CL. The input values 
(here (pt = 2 and 9t = 1) are unambiguously selected by this additional constraint. 


To recover the six unknown of the three-momenta of the two missing neutrinos one can use several 
constraints, in the narrow width(s) approximation: the two masses of the t and i and the two masses 
of the W~^ and W~. With two additional constraints one can compute the six unknown needed. In 
contrast to the pioneering work of [7|, the treatment we propose ignores possible missing-photon(s) 
and uses all four equations expressing the conservation of the four momenta. Together with the 
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four top and W mass constraints, this provides a total of eight constraints, whereas six are needed. 
The two constraints in excess are used to compute also the energies of the two 6-jets, assumed to 
be poorly measured. The directions of the two 6-jets are assumed to be reliably determined, even 
in presence of hard gluon radiation. Solving this set of eight equations with eight unknown yields 
two solutions, in a first step. These multiple solutions are further amplified when one exchanges 
the role ot 6 and 6 (assuming that the charge of the 6-jets are not determined experimentally). To 
select the right solution, one can make use of the measurement of the 6-jet energies. Although the 
6 -jet energies can be only poorly measured, compared to the lepton energies, it turns out that one 
can select the correct solution in most of the cases (see Fig. [^. 



Figure 2: The left panel shows the distribution of the reconstructed cos 9t in the ILC rest frame 
obtained with the narrow width treatment applied on Monte Carlo events produced by GRACE, 
at leading order, but taking into account the widths of the top and W. The right panel is the 
distribution of the difference between the reconstructed and true values of cos Of, in the top rest 
frame, here adjusting the top and W masses in the final fit. 


Using events generated by the GRACE [8] software, a preliminary Monte Carlo study was per¬ 
formed to start assessing the performance of the kinematical resconstruction for the fully leptonic 
case. This (leading order) Monte Carlo simulation does not account for missing-photon(s), nor for 
hadronization, nor for detector effect, but it does not rely on the narrow width approximation. 

Eor the sake of illustration, Eig. [^represents the distribution of cos 9t (the reconstructed cosine of 
the direction of t in the ILC rest frame). The same figure (right hand side) shows the distribution 
of the difference between the reconstructed and true values of cos Of, (the cosine of the direction of 
the b with respect to the top direction, in the top rest frame) using Breit-Wigner constraints to 
let the masses of top and W vary in the final stage of the fit, which improves the precision of the 
reconstruction. 

Eor about 5% of the events, the wrong 6 combination is retained. A significant fraction of the latter 
are irreducible errors: the wrong combination being indeed a better match, by accident, due to 
top’s and/or IT’s being far off-shell. 
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3 Expected sensitivity to the form factors : LO example 

3.1 Narrow Width Approximation 


In the narrow width approximation, and at born level, the distribution of events can be described 
in the helicity formalism. We denote the helicity of particle x. 


• The matrix element of the production of the tt pair through Is denoted as : 

• The matrix element of the decay t —)> b : Xt,^ 

• The matrix element of the decay i ^ l~vi b : 


For fully polarized beams, the distribution of events is described by: 

|2 


r'== E \J2m 








( 1 ) 


In the following we assume that the data analysis aims to probe the couplings of the top to the 
Z/ 7 . The coupling of the W is assumed to be the one of the Standard Model (SM): for the massless 
b, I and u only the left-handed helicities are involved, and the hrst sum can be removed. 


3.2 Form factors and the angular distribution 


By using the angular distribution of the top quark production and decays, which can be obtained 
by measuring the kinematics of the fully leptonic decays, as discussed in the previous section, we 
can extract the top quark polarization pattern. In particular, this pattern carries the information 
of the interactions which produce the top and anti-top pair. In [Hj, assuming that the tt pair is 
produced from annihilation of e'’'e“ into a spin one particle, the most general Lagrangian for the 
top and anti-top production are obtained in terms of the form factors as: 


-^int — 


v='y,Z 




( 2 ) 


Since for the EDM term, the coupling F^^ can be a complex number, we have ten real form 
factors. Note that in the literature, one can hnd different notation using the form factors F’s (see 
e.g. [in]). Our numerical result is partially given in terms of F’s for comparison but these two 
definitions are related via simple formulae: 




= -{F, 


W 

IV 


+ -^ 21 /) 


T^V _ LP'^ 

^2V ~ 


F^a = -Fia, 


F}a = -iF^A 
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By using these form factors, the angular distribution for each combination of the initial and the 
final polarization stems from the snm of the 7 and Z exchange amplitudes (see e.g. my- 


MieieR - 
MieieR - 
MicLCR 
MieiCR - 
MieRBL - 
MicRCL - 
M{eReL 
M{eReL - 


tLtRV/^ 

tRhV/^ 

■ tLhV/^ 

tLiRT'^ 

tRtLV'^ 

■ 

tR^RV^^ 


= -/3F7f+ F2>^^](l + cos0)e-*'^ 


= c 


= c 


= —c 


= c 


■'L L-^ IV lA ' 2V A" ' 

4\fIv + /5^7f + - cos0)e-*'^ 

+ ^P2a)\ sin^e"*'^ 
- /3F7f)] sin0e-*'> 
](1 - cos 0 )e*'^ 

f [F7f + /3F7f + F7f ](1 + cos0)e*'^ 
'7"'[i"7f + l\Flv + /3^2A^)] sin^e*^ 


-R 

, 7 /^ 

R 


= <^R 7 


where /3^ = 1 — Am1/s,'y = -y/s/( 2 mt) and the overall factors are: 


c7 = - 1 , 


^7 _ _i 


cf = 


- 1/2 + si 


SnnCa 


S — TTir 


Cr = 




S — nir 


(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 


( 11 ) 


where = sin 0 ^ and = cos^^y, with 6\y being the weak mixing angle. 


The core of the experimental technique advocated here is to use the complete kinematics of the 
final state to perform a likelihood analysis based on the complete PDF as provided by the matrix 
element of the process. Doing so, the analysis is optimal, since it nses all available information. 


3.3 Optimal data analysis 

Let’s define | A1 p (a) the complete matrix element squared, for a given beam polarizations. It 
depends on k constants, here collectively denoted by a. Upon integration over phase space, and 
for a given luminosity £, one expects a nnmber of events: 

N{a) = £ y I At p (a) dLips (12) 

where the Lorentz invariant phase space is: 

dLips oc dcos0i dcos^j, d(/>fe dcos^^ (lcj)i dcos0;+ d(/>;+ dcos0;- d(/;- dql dg| dg^ dg^ (13) 

where the angles are defined in the appropriate rest frames, and the g^’s are the invariant masses 
of the top’s and W’s. Note that, since the normal to the plane defined by the e"’“e“ pair and the 
ti pair can be used as a reference to measnre the azimuthal angles, all angles involved in the 6 
particle final state a priori contribute to | P, up to an overall rotation around the z axis 

that allows to set (pt = 0, which therefore does not appear in dLips. In the following illustration, 
one uses the narrow width approximation, so that the fonr g^ are integrated out. As a result the 
Lorentz invariant phase space is rednced to nine-dimension. 
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If one is using the full matrix element squared in a likelihood analysis relying only on the distri¬ 
butions of events in phase space, the inverse of the expected covariance matrix can be expressed 
as: 


where: 


V}. = Aij = N {{uji - - %))o 


N 


Ui 


fo 


N{ao) 

d\M\^ (g) 1 

dai |„o I Ai p (a*^) 


dN{a) 


1 


dai |q!oA^(q:^) 

I M p (q;°) 

/ I p (a°) dLips 


= uji fo dLips = 


(14) 

(15) 

(16) 

(17) 

(18) 


where ao stands for the set of a* values that maximize the likelihood. In the following we assume 
that they correspond to the null-hypothesis: = 0 in the SM. To take into account the constraints 

coming from the yields amounts to ignore the flj and use, instead of Eq.(14): 

= N {uiiOj 


K- 




'^ 3/0 


(19) 


The above expressions are illustrated in appendix. In the following, the analysis is assumed to 
make use of the yields (i.e. Eq.(19) is used) which contributions are important. Eor the sake of 
clarity, only the yields of the fully leptonic final state are used, and not the yields from other final 
states. 


The above Eqs.(14)-(19) convey an important set of messages: 


• For each theoretical parameter ai, and for a particular experimental set-up (e.g. for various 
beam polarizations or energies) there is a corresponding kinematical variable uJi that captures 
all the relevant information carried by each event EE]: therefore, this optimal variable should 
be accurately measured, 

• The LOi variables can be defined in any case: the matrix element does not have to depend 
linearly on the theoretical parameters, nor does it need to be cast into a readable expression 
(the LO expression is already quite intricate), 

• Different experimental set-up correspond to different Ui variables, however, to combine the 
data set is straightforward: one just needs to add the corresponding A matrix to obtain the 
full A matrix, 


It should be pointed out that the use of the Ui variables is done implicitly when one performs a 
likelihood analysis using the full matrix element squared. These variables do not bring additional 
statistical power, but they provide a means to grasp more easily what the fit is doing, and which 
events are the important ones. 
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3.4 Result using the optimal analysis method 


As a matter of illustration, if one performs a simultaneous fit of the deviations with respect to the 
Standard Model values of 10 of the form factors one gets the following error matrix (diagonal terms 
are the statistical uncertainties, while off-diagonal terms are the correlation coefficients); 


'TZe (5PjV 

7^e SFfy 

7^e 

TZe SFf^ 

TZe SF^y 

TZe SFiy 

TZe SF^^ 

TZe 6F^^ 

Im 6F^^ 

Fm SF^j^ 

0.0037 

-0.18 

-0.09 

+0.14 

+0.62 

-0.15 

0 

0 

0 

0 


0.0063 

+.14 

-0.06 

-0.13 

-t0.61 

0 

0 

0 

0 



0.0053 

-0.15 

-0.05 

-t0.09 

0 

0 

0 

0 




0.0083 

-t0.06 

-0.04 

0 

0 

0 

0 





0.0105 

-0.19 

0 

0 

0 

0 






0.0169 

0 

0 

0 

0 







0.0068 

-0.15 

0 

0 


0.0118 0 0 
0.0069 -0.17 

0.0100 

The correlation coefficients between the 6 first form factors and the last 4 form factors exactly 
vanish: this absence of correlation, comes from the property of the integral, when uii is CP even 
and ujj is CP odd (see also the appendix): 

= dLips = 0 (20) 

because f /o dLips is CP even. Similarly, the correlation between the real parts and the imaginary 
parts of and vanish because of C (and P) symmetry. 

It should be stressed that the above matrix corresponds to the ideal case; detailed studies require to 
take into account the dilution of information due to detector and physics effects, as well as related 
systematical errors. Note that, while the beam polarizations help to improve the measurements, 
they are not essential for the method to apply. 

The semi-leptonic decays, in principle, can also be analyzed using the matrix element method. In 
that case, one trades a well measured charged lepton, but a missing neutrino with two jets, not 
so well measured and for which charge information is lost (in general), and one also gains a factor 
6 in statistics. Assuming that the b — b jets can be properly assigned to the top decays, one can 
compute the matrix element squared, symmetrized with respect to the two quarks from the W 
hadronic decay. The resulting statistical errors are about a factor two smaller than for the fully 
leptonic final state, here also for a perfect detector. 



4 NLO matrix element study 


As mentioned in the introduction, the NLO corrections to the e'’'e“ —?■ tt process are large: they 
amount to ~ 5% for the total cross section and ~ 10% for the forward-backward asymmetry. Since 
the experimental measurements can reach the per-mill level, it is clear that we can actually measure 
these electroweak corrections and furthermore identify the deviation from the SM at this order. 
The GRACE program [8] can provide the SM prediction for e^e“ —tt including the full one 
loop electroweak corrections. About 150 diagrams have been computed automatically |13) . As an 
example, we present the result for the total cross section of the e'’'e“ —)• ti process in Fig. From 
the figure, we can extract the following conclusions: among large corrections of the full 0(q:) terms, 
the largest correction comes from initial state photon-radiations, which is denoted as Sqed in the 
figure. Besides the initial-state photonic-correction and another 5% of trivial weak corrections 
estimated using the scheme, still 5% corrections from non-trivial weak-correction remain. 


e e+ —)■ tt 



Center-of-mass energy [GeV] 


Figure 3: The result of the full electroweak correction and the genuine weak correction on the total cross 
section, as a function of the center-of-mass energy (taken from [l3]). The circle points represent the QED 
correction, the empty rectangle points are the results for the full electroweak correction while the triangle 
points are the results for the genuine weak correction in the a scheme. The filled rectangle points represent 
the results of the genuine weak correction in the scheme 


Recently the initial and the final state polarizations have been implemented into the GRACE 
program, which allows us to perform a more detailed study in this direction. We show two examples 
of the results obtained by this GRACE version in Fig.[^ where we investigate the NLO contributions 
for the polarized initial states, and The left panel shows the NLO contribution to the 

cross section (solid lines) for (red circle) and e~^e^ (blue square). Dashed lines are the tree 

level result and the green lines (triangle) are the sum of two polarization cases. The NLO correction 
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is significantly larger for the polarization case. The right panel shows the cos 6 dependence of 
the NLO contributions (with the same color scheme as the left panel). Interestingly, we find that in 
the case of the NLO contributions are negative (positive) for positive (negative) cos0 while 

for e^e^, the NLO contributions are always positive. This kind of strong dependence of the NLO 
terms on the kinematical variables can be most useful to investigate the NLO corrections in detail. 




Figure 4: Preliminary results produced by the GRACE software with polarized beams. The red (circle) 
and blue (square) lines are the result for initial state polarization and respectively. The solid 

lines represent the tree level plus one loop corrections while the dashed lines are tree level only. The green 
(triangle) lines are the sum of the red (circle) and blue (square) lines. The left panel is the ^/s dependence 
of the cross section and the right panel is the cos 9 dependence. Note that the drawn lines are only to guide 
your eyes. 


It is important to emphasize that the angular distribution formalism used in the previous section 
does not apply fully to the NLO electroweak corrections. For example, in Eqs.Q-(lO) it is assumed 
that e^e"*" annihilate into spin one particles, but the box diagram which appears as a one loop 
correction does not fit into this category. As a result, the angular distribution corrections from the 
box diagram cannot be written in terms of (1 ±cos 6) or sin 9 as in the above expressions. Therefore 
it is necessary to define a set of theoretical parameters that describes completely the NLO matrix 
element and use this full NLO matrix element in the analysis. Then, the numerical results obtained 
in the previous section are likely be affected. Thus, the analysis including the full NLO corrections 
is essential to be conclusive to the sensitivity to the new physics in this process. 
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5 Conclusion and future prospects 


In this paper, we showed that fully leptonic reconstruction of tt events (e+e" —)• 11 —)■ l^l~b b 
at linear colliders presents a strong potential to perform precise measurements of the tiV coupling, 
complementary to the measurements provided by the semi-leptonic process. In particular, the 
full kinematical reconstruction in this process allows us to successfully apply the matrix element 
method. We have performed a sensitivity study to ten form factors which parameterize in rather 
general terms the tiV couplings at the tree level. Based on the ILC baseline configuration, ^/s = 500 
GeV, with a luminosity of 500 fb“^, we found that the ten form factors can be simultaneously 
determined below a percent level precision. The most important next step is to include experimental 
and physics effects to obtain a more realistic assessment of the statistical power of this analysis, 
and to evaluate systematical uncertainties, both experimental and theoretical. For instance, one 
should perform the analysis using the leading order GRACE Monte Carlo producing 6-fermion final 
state, not only through ti pair production. 

We have also discussed the NLO electroweak corrections to e'’'e“ —)■ ti, known to be non-negligible 
in this process. After removing the reducible corrections (i.e. QED photon emission, corrections 
to Fermi constant), we find that there is still 5% of loop effect remaining. The high sensitivity 
demonstrated at LO by the study presented in this paper provides a proof of principle that the 
method could be applied, and encourages us to attempt an electroweak precision test using this 
process with full NLO electroweak corrections included. We have shown a preliminary result of the 
NLO computation of GRACE with polarized beams. Interestingly, we have found that the pattern 
of the NLO corrections are very different for the different configurations of the beam polarizations. 
These correlations must be important for the precise determination of the top electroweak couplings. 

Our final goal is to perform a similar analysis including fully the NLO corrections. However, a 
calculation for full one-loop corrections with six-body final states, including three-body decays of 
each top quark, is presently beyond our reach. The number of diagrams involved in this process 
is larger than 90,000 and some of them may imply more than 100GB of executable modules. 
On the other hand, the on-shell approximation for the NLO calculation is expected to ensure 
enough precision for the analysis. Eull one-loop electroweak corrections are estimated separately 
for the top-quark pair-production and for their three-body decays, then they are convoluted as 
Ale-body = ® AI(£Sy + ® . Non-factorizable corrections, such as a photon 

bridge between initial electron and final lepton, can be estimated separately, and accounted for, if 
necessary. It is important to emphasize that the form factor parameterization, thus the angular 
distributions, which we have used at the LO study would be inapplicable for the NLO study due to 
the new type of diagrams which can not be parameterized solely by the effective tiV coupling, such 
as the box diagrams. It would be very interesting to see how much the numerical results obtained 
in this article would be affected by these contributions. 
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Appendix: Illustration of the The Matrix Element Method : CP 
violation 


We consider here the special case where only a subset of CP-violating parameters are intended to 
be measured, and we derive step by step the result of Eq.(14): 


a = ar + iai = —2i sin 9w 


( 21 ) 


The matrix element can be written as the quadratic expression: 

I -MXeXe Oo + OrOir + Oitti + 02{o'l J- (22) 

where the cross term vanishes. The above expression being positive definite for any values of 
Or and (Uj, the coefficients satisfy: 

Oo > 0 (23) 

O 2 > 0 (24) 

4 O 0 O 2 > Oj + Of (25) 

In particular | O,./ i \> Oo is allowed (and does occur). Having in mind that a^/i <C 1 one may be 
tempted to ignore the contributions. However doing so does not bring any simplification, and 
it raises the concern of events having O 2 i§>| Or/i I- hence, the terms are kept in the following. 
This illustrates the fact that the use of the Ui variables is not restricted to a linear dependence of 
the matrix element on the theoretical parameters to be measured. 


Upon integration over the phase space, one gets: 

J I A^AeAe 1 ^ dLips = (Do + (DrCXr + (DiCXi J- 02(0^ T 
where the second and third term vanish: 

Or = Oi = 0 


(26) 


(27) 
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This is because Or and Oi are both odd under CP: then upon integration over the phase space, 
each point will match it CP conjugate point to cancel out the resulting integral: stated differently, 
in the case at hand the derivatives with respect to the theoretical parameters a of the yield N{a) 
vanish, for = a? = 0. 


The log-likelihood expression reads: 


L(aj) = N In 


N 

a 


In 


Oq OrOtr + OiOli O 2 {ttr “I" ^i) ^ 
Oo + 02(0:^ ) 

1 -|- OJrOir + tUiCXi -j- ^2(0:^ + of) 


- dLips 

Oo dLips + cst 


1 02(a2 + a^) 

where the constant ”cst” is irrelevant and where one introduced the notation 


Or 




UJr = 


UJi = 


Oo Oo ’ 

The inverse of the covariance matrix is given by: 


u;2 = 


O 2 

Oo 


rio — 


Oo 


Ajj — 


d‘^L{a) 
daidaj 1^^ 


Evaluated for = 0, the A^r term boils down to: 


Arr = N Ur^ dLips - ^ j O 2 dLips 202^ = N{ujI)^ 


(28) 

(29) 

(30) 

(31) 

(32) 


and similarly for the other terms. One thus recovers the result of Eq.(14), with = O^/j/Oo = 0. 
In particular only oji and ojr enter in the expression of Ajj. 


Erom the uJr and uJi distributions, as provided from the narrow width Monte Carlo simulation, one 
obtains: 

o-[a,.] = Y^A^ ~ 0.01 ; (T[ai] =-^A^ ~ 0.01 (33) 

with a limited improvement due to beam polarizations. On the other hand, we found that using 
the full matrix element, one improves the accuracy on the CP violating parameter by almost one 
order of magnitude with respect to [7], which illustrates the potential gain the method can provide. 
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